Study on the mechanical system related to Hahn's discrete time derivative 
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In this paper, we use the quantum variational calculus related to Hahn's discrete time derivative 
construct the deformed version for the classical mechanics related to the Hahn's calculus. We deal 
with the deformed dynamics such as the motion with constant velocity and the motion with constant 
acceleration. Moreover, we extend our work to the motion of a body in a resisting medium by using 
the new identity for an infinite (q, w)-series , where the retarding force is assumed to be proportional 
to the deformed average velocity. 
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q ; I. INTRODUCTION 

<D 

■ Quantum difference operators are receiving an increasing interest in the applied mathematics and the theoretical 
physics because of their application [1-5]. Briefly Speaking, a quantum calculus substitute the ordinary derivative by a 
difference operator, which allows us to treat sets of non-differentiable functions. Since Jackson [6] introduced the first 
expression of difference operator called a Jackson derivative, several expressions of the difference operator appeared. 
Among them, the most famous one is Hahn's difference operator [7], which has two deformation parameters w and q. 
Hahn's operator becomes Jackson's q-derivative when the parameter uu goes to 0. 

Recently, Malinowska and Torres [8] studied the Hahn quantum variational calculus, in which they studied the 
, Hahn quantum Lagrange problem. From now on, we will restrict our concern to the case that the parameter w goes 
to 0. Indeed, their action has the following form : 



,6 

S[x] = / C(x, D t x;t)d q t, 

J a 



(1) 



where x satisfies some boundary condition. From their quantum variation calculus, we can obtain the following 
, equation minimizing the action given in the eq.(l) : 

m : dc 



Hahn's discrete derivative is defined as 



f(qt+w)-f(t) (. / s 

D q>w f{t) = \ fa-i)t+« { l ^ w o> (3) 
' /'(*) (t = wo) 



where w = -^- q 

The Leibnitz rule for Hahn derivative is as follows 



D q , w {f{t)g{t)) = (D q . w f(t))g(t) + f(qt + w)D q . w g(t) 

Dq w (IB) = (A^/ft))ffft) - f(t)D v , w g(t) ( 4 ) 



\9(t)J g(t)g{qt + w) 

The Hahn integral is defined as the inverse operation of the Hahn derivative : 



/ f(t')d q , w t' = ((l-q)t-w)J2<l k f(t<l k +[kk W ), (5) 
Jwo k=0 
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where the (q, w)-deformed number is defined as 



[*].. = ^f 1 (6) 



The q, w-deformed number is proportional to the Jackson's g-deformed number as follows : 

[k] q , w = w[k] q , (7) 



where 



W, - *£p- (8) 



We have the following useful formula from the definition of the Halm's derivative : 



n-l 



D q , w (at + b) n = aJ2Hqt + w) + b) k (at + b) n - k - 1 (9) 



fe=0 



Inserting a = 1, b = into the eq.(9) gives 



A^ n = ^(gi + w;) fe r- fe - 1 (10) 

fc=0 

Now let us introduce the following polynomial of degree n as follows : 

n 

{t;q,w) n = Y[(t-[j] q , w ) (11) 

Then we have the following property: 

^^(i; <?, w)„ = [n]«,,o(*; 9, w)™-i (12) 

This work urges us to construct the deformed classical mechanics related to Hahn's calculus, which is the main 
purpose of this paper. 

II. (q, ui)-ANALOGUE OF GALILEI'S FORMULA 

With the formulation of the Hahn's calculus, we are able to construct the deformation of the classical mechanics. 
Let us introduce the deformed velocity v(t)&nd deformed acceleration a(t) as follows : 

v(t) = D t x(t), a(t) = D t v(t), (13) 

where x(t) is a deformed position and the deformed time derivative D t is a Hahn's derivative defined by : 

The eq.(13) is easily derived form the deformed Lagrangian defined in the eq.(l). We restrict our concern to the case 
that the w is real and q satisfies < q < 1. When w goes to and q goes to 1, a deformed time derivative reduces to 
an ordinary time derivative. 

Now let us start with the motion with the constant velocity v. The relation between the deformed position and the 
deformed velocity as follows : 



(15) 
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It can be written as 

x(qt + w)— x(t) = (q — l)vt + vw (16) 

Iterating the eq.(16) N times, we obtain 

x{q N t + [N] q w) - x(t) = {q- l)[N] q vt + [N] q wv, (17) 
where we used the following formulas : 

q k (qt + w) + [k] q w = q k+1 t + [k + l] q W (18) 

and 

N-l 

q 



Y,[k] q = ^—{N-[N] q ) (19) 



fc=0 

Inserting t — , we get 

x(w) - x(0) = vw (20) 
Now let us take the limit N — >• oo in the eq.(17). Then the eq.(17) becomes 

x{t) = x(wq) + vt — vwq (21) 

Inserting t = into the eq.(21), we get 

x(p) — x(wo) — vwo (22) 

and 

x(t)=x(0)+vt (23) 

Indeed, the solution (23) satisfies the relation (16). 

Now let us discuss the motion with the constant acceleration a. The relation between the deformed velocity and 
the constant acceleration is given by 

v(qt + w) — v(t) , . , . , H . 

'- ^=a or v(qt + w) - v(t) = (q - I) at + aw (24) 

t(q-l)+w 

It gives the following simple relation : 

v(t) = v(0) + at (25) 
Because the deformed acceleration is given by applying the Hahn's derivative to the deformed velocity, we have 

D t x(t) = v(0) + at (26) 

or 

x(qt + w)- x(t) = (t(q - 1) + w)(v(0) + at) (27) 
Iterating the relation (27) N times, we have 

JV-l 



x(q N t + [N] q w) - x(t) = J2 1 k (t(<l - 1) + w)(v(0) + a(q k t + [k] q w)) (28) 
Using the following property 

E^ = dh™^-1] 9 , (29) 
k=o 1 



k=0 



N-l 

q 



we obtain the following relation from the eq.(28) : 



x(q N t + [N] q w) - x(t) = (t(q - 1) + w)(v(0)[N] q + at[N] q 2 + -^L[N] q [N - 1],), (30) 



q 

Taking the limit TV — > oo, the cq.(30) is then given by 



= a;(«;o) - w(0)w o - + w(0) * + ^TT (31) 



q+l q 

Inserting t — into the eq.(31), we have 



2 

xtOj^W-^H-^, (32) 

q+l 

which gives the following solution : 

x(t) = x(0) + v(0)t + (33) 

III. GALILEI'S FORMULA BY USING THE SECOND ORDER (q, «j)-DIFFERENCE EQUATION. 

Because the deformed acceleration is given by applying the Hahn's derivative to the deformed position twice , we 
obtain the following relation : 

D 2 t x{t) = a (34) 
From the definition of the deformed time derivative, we can obtain the following second order difference equation : 

x(q 2 t + [2] q w) - [2} q x(qt + q) + qx(t) = qa(t(q - 1) + wf (35) 

If we set 

h(t)=x(qt + w)-qx(t), (36) 

we obtain 

h(qt + w) - h(t) = qa(t(q - 1) + w) 2 (37) 
Iterating the relation (37) N times, we get 

h(q N t + [N] q w) - h(t) = aq[N] q 2 (t(q - 1) + w) 2 (38) 
Taking the limit N — > oo, the eq.(38) is then given by 

h(t) = h(w ) + ^^(t(q-l)+w) 2 , (39) 

which gives the following : 

x(qt + w)- qx(t) = h(w ) + -^L_(i(g _ X ) + w f ( 40 ) 

This equation is still satisfied by replacing as follows : 

x(t) = x{t) + Ct - Cw , (41) 
where the constant C will be fixed by demanding the initial condition later. The eq.(41) is then given by 

x(qt + w) - qx{t) = h{w a ) + -^j(t(q - 1) + w) 2 (42) 
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Iterating the eq.(42) N times, we have 

q 1 ~ N i(q N t + [N] q w) - qx(t) = h(w ) + ^jWl ~ 1) + ™? (43) 
Inserting t = wq into the eq.(36), we obtain 

h(w ) = (1 - q)x(wo), (44) 
which enables us to replace the eq.(43) with the following equation by taking the limit TV — > oo : 

I(t, = l(0) + K' "1217' + pf < 45) 

Using x(0) — x(0) and the eq.(41), we have 

*(*) = *(«*) - C«* + ^ + (C - _ j t + —f (46) 
Inserting i = into the eq.(46), we have 

x(0) = x(wo) -Cw + ^ (47) 



Applying the Hahn's derivative to the deformed position given in the eq.(46), we have 

2 at 
"[2] 



v(t) = Aa:(t) = ( C - ) + at ( 48 ) 



Inserting t = into the eq.(48), we obtain 

r? - „fni j- : 

[2] 



C = w(0 ) + ^ (49) 



Therefore the deformed position becomes 

1 



x(t) = x(0) +v(0)t + j^at 2 , (50) 

[ 2 \q 



where 



i(0)=i(«*)-«(0H-J (51) 

IV. VERTICAL MOTION OF A BODY IN A RESISTING MEDIUM 

In the (q, w)-deformed mechanics, the (g, w)-deformed Newton's equation is easily derived from the eq.(2) as follows 

F = ma — mDtV, (52) 

where v(t) and a(t) are the q-deformed velocity and q-deformed acceleration, respectively and to is a mass of the 
body. When the force is constant, the body moves with the constant acceleration — . It should be emphasized that 
the force is not necessarily constant and, indeed, it may consist of several distinct parts. First let us discuss the one 
dimensional motion with air resistance proportional to the (g, u;)-deformed average velocity. Then the (g, w)-deformed 
equation of motion is given by 

v(t)+v(qt + w) 

mD t v = -k — , (53) 

l Z \q 
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where - »( t )+^(g t + t ") [ s called a (q, w)-deformed average velocity and becomes a instantaneous velocity when the defor- 
mation parameter q goes to 1 and w goes to 0. The parameter A: is a positive constant that specifies the strength of 
the retarding force. From the definition of the (q, w)-deformed time derivative, the eq.(53) has the following form : 

1 - -*L-(( ff - l)t + «,)) v(t) = (l + J^-((q - l)t + «,)) + w) (54) 

By iterating this functional equation N — 1 times, we obtain 

«W - VT ^«(^ + [*],«,), (55) 

where the q-shifted factorial is defined by 

(a;g) = l, (a; ?)at = (1 - a)(l - ga)(l - </ 2 a) • • • (1 - « JV " 1 a) ) n = l,2,3,--- (56) 
Iterating the functional equation infinite times gives the following : 



where we used the following relation : 



«(*) - r^v(w ), (57) 



lim = 0, lim [JVLw = w (58) 
If the initial velocity v(0) is set to vo, we have 

v(w ) = v (59) 

Then we have 

v(t) = vo^- r-2a (60) 

Now we introduce the well-known (<?, w)-deformed exponential function satisfying 

D t e q ^ w (at) = ae q>w (t), (61) 

where a is an arbitrary constant. The (q, w)-deformed exponential can be written by using the q-shifted factorial as 
follows : 

&qAat) = (-*((* -1* + «,);,)«, ^ 
Inserting the eq.(62) into the eq.(60) gives the following expression : 

v (t) = «o j— k — r- (63) 

Indeed, we can easily check that the above expression reduce to the classical result 

v(t) =v e- k ^\ (64) 

when the deformation parameter q goes to 1 and w goes to 0. 

Now we study more complicated problem. We discuss the vertical motion of a body in a resisting medium in which 
there again exists a retarding force proportional to the (q, w)-deformed average velocity. Let us consider that the 



7 



body is projected downward with zero initial velocity v(0) = in a uniform gravitational field. The (q, w)-deformed 
equation of motion is then given by 

v(t)+v(qt + w) 

mD t v = mg- k — (65) 

From the definition of the (q, u>)-deformed time derivative, the eq.(65) has the following form : 

m[2\ q J V m [ 2 J<? / 

By iterating this functional equation N — 1 times, we obtain 

v(t) = JlLvtft + [N } gW ) g((q l)t + w) £ ^ ( 67) 

The eq.(67) can be written as the following form : 



( lr((9- l)t + w);q) 

V(t) = ^V(q»t + [*],«,) 



+ 7 — — v- [ £ ] ? , (g ; g)A ; — 9 " (2n+i) f 4r«ff - 1)* + *>)) 2n+1 > ( 6s ) 



(i((«- 1 )* + «');?) JV ™=° 



where [a;] is a Gauss symbol. Because the initial velocity u(0) is zero and limAr^oo q N = 0, we have following relation 
by taking the limit N — > oo : 



«(*) = -7—7 H^K) 



hm ? PW/* ^ ± i^-f^^ f4r(-(t- W Q))) 2 " + \ (69) 

JV_KX, (i((«- 1 )* + «');«) JV (?;g)2„ + i(?;g)iv-i-2„ \m[2] g 



Because (q;q)N and (q; g)jv-i-2n reduce to 1 when TV goes to the infinity, the expression (32) takes the following 
form : 



v{t) = ^H»o) 

+ 7 [2] ^ g/fc N f (i^v^ ( * ( _ (t _ ^ 2 " +1 m 



m[2], 



or 



= t ^f(wo) 



3.™ (mp],*) 



[2\qmg 
k q ' w 



m[2] q 



'I 



(1 - g)2n+l g n(2n+l) / fc 



n=0 



(9; 9) 



2n+l 



■(-(*- «*>)) 



2n+l 



The eq.(71) can be written by using the q-number factorial as follows : 



v(t) = 



v(w ) + 



[2] q mg 



k eq,w 



^ _n(2n+l) / fc 



™[2], /^[2n + l],! Vm[2] 



-(-(*-«*>)) 



2n+ 



where the q-number factorial is defined by 

[n] q \ = [n] q [n-l] q ---[2] q [l] q 
Comparing the eq.(72) with the definition of the q-deformed exponential function such as 



00 1 



we can rewrite the eq.(72) as 



v{t) 



<1> W \ m[2] q 



where we used the following relation 



q,w 



E 



-(-(*- t«o)) 



1 - q~ n 

[n\ q -i\ = q 2 [ n ] ? ! an d [n\ q -i = - - 
Thus we can express the eq.(76) by using the q-deformed exponential function as follows : 



«(*) = 



e 9>«> ( mp],*) 



,(-0) + fee D 



where we used the following identity : 



2k 



m[2] q 



m[2] q 



t - e„ 



m[2] q 



-,2n+l 



e 9 - 1 W-e 9 - 1 (-a) = 2g I ^ 1]9 _ i! 



Inserting t = into the eq.(77), we have 
Then the equation (77) is written by 



f (0) = v(w ) 



v(t) = v(0)- 



2fc ' 



»[2] s 



t - e„-i 



'[2],' 



If the initial velocity is set to 0, we have 



m[2] q 



m[2] q 



2k ~ q ' w 

It can be easily checked that the expression (77) reduce to 

t,(t) = ^(l-e-«M), 

when he deformation parameter q goes to 1. 



t - e„-i 



m[2] q 
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V. CONCLUSION 

In this paper, we used the quantum variational calculus related to the Hahn's calculus [8] to construct the (q, w)- 
deformed version for the classical mechanics. We dealt with the motion with constant velocity and the motion with 
constant acceleration by using the first order (q, w)- difference equation. We also obtained the same result by using 
the second order (q, w)- difference equation Moreover, we extended our work to the motion of a body in a resisting 
medium, where the retarding force is assumed to be proportional to the (q, w)-deformed average velocity. In this 
process, we found the new identity for infinite q-series. 

We think that our work will be extended to more complicated mechanical problem such as (q, u;)-deformed damped 
oscillator problem , (q, w)-deformed Kepler problem and so on. We hope that these topics will be solved in the near 
future. 
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